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ABSTRACT 

We combine radial velocities, proper motions, and low resolution abundances for a sample of 315 K 
and M giants in the Baade's Window (I, b) = (0.9°, —4°) Galactic bulge field. The velocity ellipsoid of 
stars with [Fe/H] > —0.5 dex shows a vertex deviation in the plot of radial versus transverse velocity, 
consistent with that expected from a population with orbits supporting a bar. We demonstrate that 
the significance of this vertex deviation using non-parametric rank correlation statistic is > 99%. 
The velocity ellipsoid for the metal poor ([FeH] < —0.5) part of the population shows no vertex 
deviation and is consistent with an isotropic, oblate rotating population. We find no evidence for 
kinematic subgroups, but there is a mild tendency for the vertical velocity dispersion erf, to decrease 
with increasing metallicity. 

Subject headings: Galaxy: bulge — Galaxy: kinematics and dynamics - Galaxy: abundances 
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1. INTRODUCTION 

The galactic bulge is the nearest example of a 
bulge/spheroidal population. A factor of 100 closer than 
M31, the proximity of the bulge not only permits de- 
tailed measurement of stellar abundances and radial ve- 
locities, but also of stellar proper motions. The first such 
study undertaken (at th e suggestion RMR) was that of 
ISpaenhauer et all (1992). The proper motions were mea- 
sured from plates obtained at the 200-inch prime focus in 
1950 and 1983; these data produced measurable proper 
motio ns on the order of .5 ma syr -1 . 

The fSoacnha uer et all (|1992f ) study was the source of 
an input catalog for a subsequent study of bulge radial 
veloc i ty and low resoluti on abundances (|Terndrup et al.l 
Il995t ISadler et al.l Il996h . An early study of the abun- 
dances and kine matics of th ese stars, using published 
abundances from iRichl (|1988| ) and radial v e locitie s from 
I Rich! (|1990l ) was undertaken bv lZhao et ail (|1994 ). That 
study suffered from having a relatively small sample size, 
and the lack of a definitive bulge iron abundance scale; 
nonetheless the most metal rich subset of stars exhib- 
ited a vertex deviation of the velocity ellipsoid in the 
plot of rad i al ver sus transverse velocity. Modeling by 
I Zhao et all (|1994l ) showed that such a vertex deviation, 
or rotation of the major axis of the velocity ellipsoid, is 
a feature expected from the long axis orbits that would 
su pport a bar stru c ture. Such a bar had been claimed 
bv lBlitz fc Spergel (|1991f ) based on asym metries in sur- 
face p hotometry of the bulge region, and bv lBinnev et al.l 
(|1991[ ) from the kinematics of molecular gas. 

The possibility of a dominant stellar bar population 



was strengthene d by the spectacu lar results from the 
COBE satellite (|Dwek et all 11994) . All-sky maps of 
the Milky Way showed an asymmetry that was con- 
sistent with a bar, with major axis in the first quan- 
drant. Many subsequent studies confirm the presence of 
a bar structure based on deproj ecting the spatial d istri- 
bution of the stellar population ( Bin nev et al.lll997l ) and 
the Schwarzschild method has been used to build a self 
consi stent rapidly rotating bar {Zhao 1996: Hafne r~et al.l 
120001 

Here we co mbi ne the proper motions of 
I Spaenhauer et al.l (11992ft . the abundances of 
i Terndrup et al.l (I1995D . and the radial velocities of 
ISadler et al.l ( 19961 ) to explore the correlation of abun- 
dances and kinematics i n the bulge stell ar population. 
Surveys now in progress (|Sumi et al.ll2004f ) have achieved 
total numbers of order 5 x 10 6 stars. However, the 
sample of 315 stars on which we report here is at present 
the largest that has proper motions, radial velocities, 
and spectroscopically estimated abundances. 

In exploring correlations of abundances and kinemat- 
ics, our aim is to constrain models for the origin of the 
bulge and bar. Wh en t he composition of bulg e stars 
(jMcWilliam fc Richlll994 : iFulbright et"aT1l2006b( ) is con- 
sidered, it c an be argued that the bulge formed earl y 
and rapidly dMatteucci et al.lll999l : iFerreras et al.ll2003D . 
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iKuiiken fc Rich! ( 20021 ) iso lated the old bulge popula- 
tion by proper motion, and IZoccali et al.l (|2003f ) demon- 
strated a comparison with a 13 Gyr old main sequence 
turnoff population; these studies find little room for a de- 
tectable trace young population in the bulge. However, 
there are complicating issues in the chemical evolution. 
Mg is found t o be elevated in bulge stars, while O is 
only mildly so (|Mc William fc Rich|[l99l IFulbright et all 
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I2006U iLecureur et alj |2006) and the "explosive alphas" 
Ca, Si, and Ti are also less enhanced than Mg, tracking 
the thick disk trends with [Fe/H]. Finally, the bar poses 
a problem for early formation and enrichment, as the 
longevity of bars is an issue. If the metals in the bulge 
were built early and rapidly, how and when did the bar 
appear and why has it survived? 

Correlating abundances and kinematics gives another 
viewpoint on chemical evolution. Even if the abundance 
distribution appears to fit the closed box model and the 
population is old, the correlation of abundances and kine- 
matics may reveal distinct kinematic subgroups. Adding 
two proper motion vectors to the radial velocity increases 
the power of this analysis. There have been indica- 
tions of abu ndance/kinematics correlations in the bulge. 
iRichl (|1990f ) found that the m ost metal rich stars have 
a smaller velocity dispersion. ISpaenhauer et al.1 (119921) 
found a similar trend in their small sample. I Zhao et all 
(1994) found evidence for vertex deviation in the most 
metal rich stars. In a rad ial velocity study of a large 
number of bulge M giants, ISharples et al.1 (1990) found 
indications for a smaller velocity disper sion for a metal 
rich subgroup of stars. iMinnitil (|1996l ) obtained spec- 
troscopy in a bulge field at (l,b)=(8,7) and found differing 
kinematics (more rotation support and a colder velocity 
dispersion) for stars with [Fe/H] > —1. In this paper, we 
explore abundance/kinematics correlations using both 
the radial velocity and proper motion data, and abun- 
dances from low resolut ion spectroscopy . We us e the new 
iron abun dance sca le of Fulbright et al. (2006a|) to recal- 
ibrate the I Sadler et alj (1996) low resolution abundances 
(those being based on the average of the Fe 5270 and 
5335 features in the spectra). While the derived abun- 
dance for any individual star is uncertain to roughly +0.3 
dex, the calibration ensures that the kinematic transi- 
tions we report here are tied to the best currently avail- 
able iron abundanc e scale. The most metal rich stars in 
I Sadler et~a l. ( 1996) have formal metallicities in excess of 
1 dex; we do not take those values seriously, as they are 
transformed from outliers. 

In this Letter we describe our kinematic analysis of the 
sample of 315 red giants. 

2. ANALYSIS 

In an axisymmetric potential, one of the axes of the 
velocity dispersion tensor of an equilibrium population 
must lie in the meridional plane. Any deviation from this 
situation indicates a non-axisymmetry, or a time varia- 
tion. We will use the space velocities of our stars to mea- 
sure the orientation of the axis of the velocity dispersion 
tensor in the radial-longitudinal velocity plane (known 
as the vertex deviation l v ) and determine whether it is 
consistent with zero. The angle l v is defined as 

tan(2Z„) = 2 ^ 2 (1) 

\ a r a l I 

where cr; r , ay and ai are the covariance and standard de- 
viation of the measured velocity components along the r 
(line of sight) and I (galactic longitude) directions. We 
measure the dispersions with an iterative clipping algo- 
rithm, by repeatedly rejecting stars that lie outside the 
ellipse with principal axes of 2.5 times the correspond- 
ing measured dispersion. Table [T] and Fig. [T] summarize 
the velocity ellipsoid parameters found, both for the full 



TABLE 1 
Velocity ellipsoid parameters.. 



[Fe/flj 


ay cr ; 
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lv 


r s 


N (N roj ) 


Prob rg 




km/sec 




deg 
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All 


109 94 


86 


-25 


-0.15 


315(21) 


0.85 


<-0.5 


108 85 


94 


28 


+0.09 


77(11) 


39 


>-0.5 


107 96 


83 


-30 


-0.23 


238(11) 


0.02 



sample and for metal poor and metal rich subsets. For 
each subsample we have calculated the Spearman rank 
correlation coefficient r s and its significance. 

The most significant vertex deviation is seen in Fig. [T] 
for the metal rich sample ([Fe/H] > —0.5), whereas metal 
poor stars show no significant deviation. This r esult is 
consistent with the results from IZhao et all (|1994f ). based 
on ~ 20% of the present sample, and proves the re- 
liability of their results (Zhao et al. also saw a vertex 
deviation in the metal-poor sample, of opposite sign to 
the metal-rich set: with the addition of the new stars in 
this study the significance of this deviation has dropped) . 
This difference between metal rich and metal poor stars 
is a strong stellar-dynamical demonstratio n of the exis- 
tence of a stellar bar in the galactic center (|Binnev et al.l 
119911 : IZhao et a l. 1994), and suggests that much of the 
disk takes part in it. 

Fig. [2] shows the velocity dispersion components and 
vertex deviation as a function of the metallicity, using 
a running mean over metallicity bins of width 0.4dex. 
The Figure clearly shows a rather sharp transition near 
[Fe/H] = —0.5: more metal-poor stars show little or no 
deviation, whereas the more metal-rich stars' velocity el- 
lipsoid is misaligned by about 30°, with a rank correla- 
tion between radial and longitudinal velocity components 
for these stars that is significant at the 99% level. A fur- 
ther indication of this change in kinematic properties is 
the dispersion in the (completely independent) vertical 
velocity component, ijb, which shows a clear break at 
[Fe/H] = —0.5 with the more metal-poor stars' b veloc- 
ity dispersion being higher. 

We have also examined the correlations between the b 
velocity and the other components, but found no signifi- 
cant deviation of the velocity ellipsoids out of the plane. 

The giant stars thus consist of a metal-rich population 
with a significant vertex deviation of 30 degrees, and a 
metal-poor population that has no significant deviation, 
and hotter kinematics. 

The general trend of the bulge to have this vertex de- 
viation is also seen in three minor axis fields with proper 
motions from HST imaging with WFPC2 (Soto et al. 
2007 in prep.). These fields were observed for radial ve- 
locities using the integral field spectrograph on VIMOS 
at VLT. The combination of radial velocities and I proper 
motions yields an observable vertex deviation in the fields 
(in samples not segregated by metallicity). 

3. DISCUSSION 

Our observations support the presence of a striking 
transition in the kinematics of the bulge, from an ap- 
parently isotropic oblate spheroid to a bar, at [Fe/H] = 
—0.5. We have used a non-parametric rank correlation 
statistic to determine that the significance of the vertex 
deviation for the metal-rich stars is > 99%. There is a 
corresponding but shallower trend toward lower vertical 
velocity dispersion (ab) at higher abundance. It is dif- 
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Fig. 1. — Velocity ellipsoids for metal rich and metal poor populations in Baade's Window. An iterative clipping algorithm to reject stars 
outside the velocity ellipsoid at 2.5 cr has been applied. Left to right, full sample, stars with [Fc/H] < —0.5, and stars with [Fe/H] > —0.5. 
In each panel the vertex deviation l v , the Spearman rank correlation coefficient rg, and the number of stars N, is listed. 
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Fig. 2. — The kinematics as a function of metallicity. Top row: Left, Vertex angle l v vs. [Fe/H] using a running box car of width 0.4 
dex. Center, the corresponding significance of the vertex deviation, using a non-parametric Spearman rank correlation. The combined 
significance for all stars with [Fe/H] > —0.5 is greater than 99%. Right, histogram showing the number of stars in 0.2 dex [Fe/H] bins. 
Bottom row, velocity dispersions in b, I and r, for the same running metallicity bins. Note the break in the kinematics in trt, at [Fe/H] = —0.5. 



ficult to imagine how a system with such a kinematic 
transition could be considered to have evolved in a one- 
zone chemical evolution scenario. 

Do abundance trends of any other elements show a 
tr ansition at [Fe/H] = —0 .5? Examining the trends g iven 
in lFulbright et all (|2006bf ) and lLecureur et alj (|2006f i. we 
find no break in the trends of Mg, O, Ca, Si, Ti, Na, and 
[Al/Fe] with [Fe/H] tha t would correspond to —0.5 dex. 
iFulbright et all (|200 6b) consider the trend of the explo- 
sive alpha elements (Ca + Si + Ti) = a cxp with [Fe/H]. 



They find that even to the lower limit of bulge [Fe/H], 
at [Fe/H] < —0.5, the bulge trend of [a ea; p/Fe] lies above 
the locus of halo and thick disk stars. Hence, the stars 
with [Fc/H] < —0.5 cannot be considered simply to be- 
long to the halo or thick disk. Similarly, virtually all 
stars with [Fe/H] > —0.5 have strongly enhanced Mg, a 
strong indication that the chemistry reflects a history of 
rapid enrichment that extends through the boundary of 
this kinematic transition. 
In order to form the bar via a stellar-dynamical in- 
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stability, which then thickens into a bulge (|Raha et al.l 
Il991| ). the stars would already have to have been in place. 
The kinematics thus suggest a scenario in which all stars 
with [Fe/H] > —0.5 stars were formed in a disk before the 
bar instability occurred; the more metal-poor stars are 
then a population unrelated to the disk, bar and bulge, 
consistent with their hotter kinematics. But the con- 
tinuity in bulge abundance trends across [Fe/H] = — 0.5 
argues against the existence of two entirely separate pop- 
ulations with different chemical evolution histories. Ev- 
idently the situation is more complicated than a simple 



few-component model can describe. It is clearly impor- 
tant to explore the existence of these trends in larger 
samples and in additional fields in the bulge. 
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